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Abstract 

It is shown that the topological discrete sine-Gordon system introduced by Speight 
and Ward models the dynamics of an infinite uniform chain of electric dipoles constrained 
to rotate in a plane containing the chain. Such a chain admits a novel type of static kink 
solution which may occupy any position relative to the spatial lattice and experiences no 
Peierls-Nabarro barrier. Consequently the dynamics of a single kink is highly continuum 
like, despite the strongly discrete nature of the model. Static multikinks and kink- 
antikink pairs are constructed, and it is shown that all such static solutions are unstable. 
Exact propagating kinks are sought numerically using the pseudo-spectral method, but 
it is found that none exist, except, perhaps, at very low speed. 

PACS Classification: 05.45.-a, 05.45.Yv, 05.50.+q 

1 Introduction 

Nonlinear Klein-Gordon kinks are found in many branches of theoretical physics. For appli- 
cations in condensed matter and biophysics, the nonlinear Klein-Gordon equation is usually 
a continuum approximation of a more fundamental spatially discrete system: kinks really 
live on a lattice reflecting some crystalline or biomolecular structure. It has long been recog- 
nized that spatial discreteness introduces important new effects into the dynamics of kinks 1J. 
Most notable is the Peierls-Nabarro (PN) barrier. Due to the loss of continuous translation 
symmetry, static kinks may generically be centred only on a lattice site or exactly half way 
between lattice sites. The PN barrier is the energy difference between these two solutions. 
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A kink moving with too little kinetic energy cannot surmount the PN barrier and becomes 
trapped. A kink moving with high kinetic energy loses energy by emitting radiation and may 
also, eventually, become trapped. So the dynamics of discrete kinks is much more complicated 
than their continuum counterparts. 

This rich behaviour has led mathematical physicists to study discrete Klein-Gordon sys- 
tems extensively for their own sake, independent of any physical application. In this regard, 
there has been a recent resurgence of interest in so-called exceptional discretizations of non- 
linear Klein-Gordon models: those which, despite the loss of translation symmetry, maintain 
a continuous family of static kink solutions, centred anywhere relative to the lattice (loosely, 
a continuous translation orbit of static kinks). An early example of this was the topological 
discrete sine-Gordon (TDSG) system j2j, which eliminated the PN barrier by preserving the 
Bogomol'nyi bound on kink energy. The idea was subsequently applied to the 4 model [3| and 
then generalized to all Klein-Gordon models [I]. A rather different approach is to choose the 
continuous family of static kinks in advance, then reverse-engineer a discrete Klein-Gordon 
system which supports them the so-called inverse method. Both this method and the 
topological discretization approach discretize the Lagrangian of the system, so they preserve 
an energy conservation law. Kevrekidis has given a method for constructing exceptional dis- 
cretizations of the PDE itself, which destroys the variational set-up, but preserves a momentum 
conservation law |HJ. With collaborators, he has recently shown that, if one discretizes the 
Laplacian in standard fashion, it is impossible to preserve both energy and momentum con- 
servation Note that in the absence of a conserved energy it is not meaningful to speak of 
a PN potential. Most recently, Barashenkov, Oxtoby and Pelinovski [S] have systematically 
constructed all polynomial exceptional discretizations of the 4 PDE with standard discrete 
Laplacian, finding families which generalize the discrete 4 systems of [3| and jSj. The new 
exceptional discretizations they find preserve neither energy nor momentum conservation. 
Moving outside the polynomial class, Dmitriev, Kevrekidis and Yoshikawa have constructed 
yet more exceptional discrete Klein-Gordon models, concentrating on the 4 case, some of 
which conserve energy, some momentum, but none both 

All these systems are theoretically interesting, but from a physical viewpoint they look 
rather contrived, and it is hard to imagine a concrete physical system that they model. The 
purpose of this paper is to point out that one of them does, in fact, model a simple physical 
system: a uniform chain of electric dipoles. The kinks interpolate between the two vacuum 
states where the dipoles are aligned uniformly forward and uniformly backward along the 
chain. This system has three parameters: the dipole strength d, the inter-dipole separation 
S and the dipole moment of inertia I. For all values of these parameters, this system can 
be reduced to the TDSG system at a single fixed value of the lattice spacing (h = VT2, it 
turns out). This allows a concrete physical reinterpretation of much of the previous work 
on the TDSG system GUI El E] It also motivates us to investigate the TDSG system 
more deeply, yielding some new results. We will construct static multikink and kink-antikink 
pair solutions, and prove that all static solutions apart from the kink, antikink and vacua are 
unstable. We will also numerically seek exact propagating kinks using the pseudo-spectral 
method, but will find only nanopterons: travelling kinks with spatially oscillatory tails. The 
tail amplitude becomes extremely small for low propagation speed, so that the possibility of 
exact propagating low speed kinks cannot yet be discounted. 
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In section we will introduce the dipole chain and demonstrate that it is modelled by 
the TDSG system. In section |3] we will construct static solutions with arbitrary winding 
number and prove that all "non-Bogomol'nyi" static solutions (that is, all solutions except 
kinks, antikinks and vacua) are unstable. Section 0] presents our results on exact propagating 
solutions, while section contains some concluding remarks. 



2 Dipole chains 



If two electric dipoles of moments d 1; d 2 are held at positions ri, r 2 G M 3 , the potential energy 
due to their electrostatic interaction is 

XT? 1 di-d 2 -3(n-di)(n-d 2 ) , , 

° E =~ A i 13 > I 2 - 1 ) 

Atteo |ri - r 2 | d 

where n = (ri — r 2 )/|ri — r 2 |. Consider an infinite chain of equally spaced dipoles, all of equal 
strength, each free to rotate in some fixed plane containing the axis of the chain. Without 
loss of generality we can assume the chain is directed along the x-axis and the dipoles rotate 
in the (x, y)-plane. Then the nth dipole is located at r n = (5n, 0,0) and has moment d n = 
d(cos0 n , sin0 n , 0), where 5 is the inter-dipole spacing and cj) n is the angle between d n and the 
x-axis. Hence, the potential energy of a neighbouring pair d n , d n+1 is 

d 2 

SE n = -^rjp [cos(0 n+ i - <f> n ) - 3 cos 4> n cos 4> n+1 ] 

d 2 r i i 

= ± p Sin2 2 ^ n+1 - + 3 Sin2 g ^ n+1 +( f>n)- 



(2.2) 



It is convenient to add d 2 /2keq5 3 to SE n so that the minimum energy of a neighbouring pair 
is normalized to 0, occurring when d n = d n+1 = (±<i, 0,0). The total potential energy of the 
chain, neglecting longer range interactions, is 

d 2 { 1 1 1 

E P = 3 ) sin2 2 ^ n+1 ~~ ^ + 3 sin2 2 ^ n+1 + f ' ( 2 - 3 ) 

while its total kinetic energy is 

^sH^y. <2 - 4) 

where / is the moment of inertia of a dipole. The Lagrangian governing its time evolution is 
then L = Ek — Ep. 

If we define the natural units of energy and time to be E = \/3d 2 /27r£ <5 3 and T = 
(I / \^3E )i respectively, this Lagrangian coincides with that of the topological discrete sine- 
Gordon system |2j with lattice spacing yl2. More precisely, let us define rescaled energy and 
time variables Ek = Ek/Eq, Ep = Ep/ Eq and t = t/To- Then 

E K -E P ='iy]{<f? n - (^^~{<f>r*i-<f> n )\ - sin 2 i(0n+i + <f>n) } (2.5) 
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where h = yl2 and / := df/dt. This is precisely the TDSG Lagrangian. It is convenient to 
take S as the unit of length. Note that 5, d, I affect the energy and length scales of the system, 
but not its dynamical properties. These depend only on h, which is fixed at y/l2 for all dipole 
chains. Note also that the physical lattice spacing of the chain is 5, not h. We shall henceforth 
use the rescaled variables exclusively, and drop the ~ superscripts. 
The equation of motion of the TDSG system is 

4 - h 2 4 + h 2 

cos0 n (sin0 n+ i + sin0 n _i) — — sin0 n (cos0 n+ i + cos0 n _i). (2.6) 



Ah 2 

If < h < 2, it supports static kink solutions 

K = 2tan" 1 e^~ b \ a = log , (2-7) 

centred anywhere on the lattice (the kink position b may take any real value). These kinks 
all have energy 1, the minimum possible for a configuration with kink boundary behaviour 
(lim^ioo <f) n = k±7r with \k+ — fc_| = 1). Thus they experience no PN barrier - this is an 
exceptional discretization of the continuum sine-Gordon model, which we recover in the limit 
h — > 0. Note that adding any integer multiple of 2tt to a single 4> n preserves solutions of (|2.6)L 
so one should think of <p n as an angular coordinate of period 27r. The constant sequences n = 
and <p n = it are two distinct vacuum solutions, between which the kink (J2.7j) interpolates. The 
dynamics of a single kink was studied in detail in It was found that kinks may propagate 
with arbitrarily low speed and never get pinned. Fast kinks on coarse lattices (h close to 2) 
do excite significant radiation in their wake which causes them to decelerate appreciably, but 
this effect is far weaker than in the conventional discrete system with the same h. We shall 
return to the question of whether the system admits exact propagating kinks with constant 
speed in section |3J but in any case kinks are very highly mobile in this system. Kink-antikink 
collisions were studied numerically by Kotecha in JU] > i n the case where the kink and antikink 
are fired at one another with equal speed. Such collisions were found to be rather similar to 
continuum 4 kink-antikink collisions. There is a critical speed v c , depending on h, above 
which the kink and antikink pass through one another. Kotecha found that v c , measured in 
lattice cells traversed per unit time, depends approximately linearly on h. Below v c there is 
a complicated (possibly fractal) structure of velocity windows in which the kink and antikink 
"bounce" off one another B > 1 times, where a bounce consists of the kink passing through 
the antikink, then turning around and passing back through in the opposite direction. There 
seem to be windows with B = oo, meaning that the kink-antikink pair settles into a breather- 
like oscillatory bound state. In fact, the TDSG system supports exact breather solutions [TTj . 
at least for h close to 2, though these are all unstable [T2*] . 

In the case of the dipole chain, h = \[V1 > 2, so the results of (3110111111121 do not directly 
apply. To make use of them we must exploit a symmetry of f)2.fij) : <f> n (t) is a solution of the 
spacing h system if and only if 

1> n (t) = (-l)> n (2t//i) (2.8) 

is a solution of the spacing h* = 4/h system. Note that if h > 2 then < h* < 2. We can thus 
identify the dipole chain with the TDSG system with lattice spacing h* = A/y/l2 = 2/y/3 < 2, 
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so the dipole chain supports a continuous translation orbit of "alternating" kinks 



<p n = (-l^tan-VM a = log [y^j > ( 2 - 9 ) 

parametrized by 6 G M, all of which have energy 1. Note that <p n still has kink boundary 
behaviour because <p n is angular, so n and —n are identified. The kink experiences no PN 
barrier. These kinks are very highly discrete, in that their structure is spread over very few 
lattice sites. For example, 86.6% of the energy of the 6 = kink resides in the — 1 < n < 1 
section of the chain. This kink is depicted in figure H 




-4-3-2-1 o 1 2 3 



Figure 1: A site centred dipole kink (6 = 0). 

Since h* rs 1.15 is close to 1, the numerical results of show that its propagation at 
modest speed is strongly continuum like: radiative deceleration is weak and the kink speed 
varies very little over a single lattice cell transit. Like their h < 2 counterparts, therefore, 
kinks in dipole chains are highly mobile. Radiative deceleration is essentially negligible if the 
frequency with which the kink centre hits lattice sites is below the bottom edge of the phonon 
band of the lattice [2]. In the case of dipole kinks, this happens when 

2 1 1 , , , s 

2ttv < - = — ^ v < -= f» 0.092 sites/unit time. (2.10) 

h y3 27rv3 
In "laboratory" units, this becomes 

v< - d (2.11) 

Turning to kink-antikink collisions, we can estimate the critical incidence speed, above which 
they pass through one another, for the spacing = 2/y/3 lattice, using Kotecha's linear fit: 
v c ~ 0.120 sites/unit time ^U]- From (J2.8)) it then follows that the critical speed for dipole 
kink-antikink collisions is 

v c w 0.208 sites/unit time w 0.144 j (2.12) 

Ve I5 

The numerical results of jTTj show that is too far from 2 for the dipole chain to support 
any breather solution constructed by continuation from the anti-continuum limit. 
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3 Static solutions 



An important consequence of the PN potential in conventional discrete systems is that they 
can, in contrast to their continuum counterparts, support finite energy static multikink solu- 
tions. The point is that if two kinks are well separated, their mutual repulsion is too weak 
to overcome the pinning force due to the PN barrier. One can thus construct static solutions 
with arbitrarily many kinks and antikinks in stable equilibrium. Since the TDSG system has 
no PN barrier one might expect it to have no static multikinks either. However, this turns 
out to be false: one can construct static solutions with lim^-too n = k±ir for any integers 
k + , k_. Only the vacua and (anti)kink are stable, however. We shall assume < h < 2, so 
the results apply to the dipole chain if we choose h = 2 /a/3 and apply ([2.8)1 . 

It is helpful to think of the static field equation as the variational equation for the functional 
E P . To be precise, for any k_,k + G Z, let 



i 2 k _ M := {0 : Z - R | J2 (*» - k -*? + X>» ~ < °°>" ^ 

n=— 1 n=0 

Note that £q = £ 2 , the Hilbert space of square summable sequences, and £k_,k+ is an affine 
space modelled on £ 2 , whose elements have boundary behaviour lim n ^± OQ <f) n = k±7i. It is 
straightforward to verify that Ep : £\_ k+ — > R is C 2 , and that is a critical point of Ep if 
and only if each of its constituent triples (0 n -i> 4>m 4>n+i) satisfies the static field equation 

dE P A-h 2 . . A — )— h 2 , , . . . 

— — = — — cos n (sm n+1 + sin <p n _{) H — — sm n (cos (p n+1 + cos <p n _{) = 0. (3.2) 

u(p n on on 

So critical points of Ep on l\_ k+ are static solutions of the model with lim n ^ = |- 00 = k±ir, and 
vice versa. Note that any triple with n +i — <j) n -i = vrmod27r automatically satisfies ()3.2|1 
irrespective of its central value <p n , by 7r-antiperiodicity of cos and sin. Such triples turn out 
to be important for our analysis of static solutions, so we make the following definition: 

Definition 1 A consecutive triple (0 n -i> 0n> <?Wi) i n a sequence <fi is called exceptional if 
4> n +i — <Pn-i — (2k + I) 71 " f° r some k G Z. 

To discuss stability of static solutions, we must introduce the second variation, or Hessian, 
of Ep. Let be a critical point of Ep and s '* be a two-parameter variation of it, that is, a 
smooth map (— e, e) x (— e, e) — > i%_ k with 0,0 = 0. Let u = <9 s s '*|( O) o), v — dt4> s,t \(o,o) ^ 
Then the Hessian of Ep at is the symmetric bilinear form Hess<i : £ 2 x £ 2 — > R defined by 



d 2 

Hess^(M, v) 



dsdt 



Ep[<f> a '*\. (3.3) 

s=t=0 



We say that is stable if the associated quadratic form u \— > Hess^UjU) is positive semi- 
definite (i.e. Hess^(M, u) > for all u), and unstable otherwise. This definition is motivated by 
the analogy with a point particle moving in M 3 under the influence of a potential V, since then 
saddle points and maxima of V are clearly unstable equilibria. We will consider the linear 
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stability criterion later in this section. By considering the variations s '* = + (s + t)u and 

0* = + tu, one sees that 

u f \ # E pW 

Hess^(M, u) = 



dt 2 



(3.4) 

4=0 



so to prove instability it suffices to exhibit a one-parameter variation 0* of such that 
(PEp^/dt 2 < at t = 0. 

Let us briefly recall the construction of the kink solutions of the model (J2.7|) . The key idea 
is to write Ep in the form 

Ep[4>\ = \ J2( D n + A* := ~ sin ^(0 n+1 - n ), F n := sin ~(0 n+1 + n ), (3.5) 

n 

and observe that -D n F n = — Acos0 n , where Af n := (f n +i — f n )/h. Hence 

< ^ ^(£>n - F n ) 2 = E P + ^ ^ A cos 4> n = E P + ^(cos k+TT - cos A;_tt). (3.6) 

n n 

It follows that Ep > 1 on £qi> an d -^p = 1 if and only if Z} n = F n for all n. This discrete 
Bogomol'nyi equation may be written 

n+ i 2 + h 6 n 

tan = -tan — , (3.7) 

2 2- /i 2 ' v ; 

and its general non- vacuum solution is (|2.7)1 . On £ 2 , a similar argument beginning with 
summand (D n + F n ) 2 shows that E P > 1 and E P = 1 if and only if _D n = — F n for all n, that 
is, 

n+ i 2-h 6 n 

tan = -tan — , (3.8) 

2 2 + /i 2 v ; 

whose general non- vacuum solution is 

n = 2tan- 1 e- a(n - b) . (3.9) 

We shall refer to (|2.7)L ()3.9|) and the vacua as Bogomol'nyi solutions of the TDSG system. 
Since they globally minimize Ep within their boundary class, they are automatically critical 
points of Ep with non-negative Hessian, and hence stable static solutions of the model. 

To construct non- Bogomol'nyi static solutions, we exploit the possibility that may have 
exceptional triples. An interesting fact about such a triple is that the potential energy con- 
tributed by its pair of constituent links is independent of the values of the triple. 



4+h 
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Lemma 2 The potential energy contributed by any exceptional triple (0 n -i, 0™, 0n+i) is Ui 
Proof: Since n +i = n -i + (2A; + l)vr for some k G Z, the energy of the pair of links 

(0 n _l,0„), (0„,0 n+ i) iS 

^(Dl-l+Dl + Fn-X+Fn) = ^ \ " ^n-l) + Sm 2 ± (0 n _ X + 7T - 0„)] 

ft, r „ 1, . 2 1 . 4 + /l 2 

+ ^l sm ^{(Pn + (Pn-l)+Sm ~(0n-i + 7T + n )j = 
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□ 



We now have 

Proposition 3 Let e C,\_ k+ be a static solution with energy E and k be an integer. Then 
<f> G ^_,fc++i de fi n ed by 

h n < k 

K-2 + 71" n> k + 1 

is also a static solution, with energy E = E + . 

Proof: Since is a static solution, every triple (0 n -i, 0n, 4>n+i) with n =<£ {k, k+1} satisfies 
(13. 2j) . while fc+1 — = fc+2 — 0fc = tt by construction, so the triples centred on n = k and 
n = k + 1 are exceptional. Hence is a static solution. 

To compute the energy of note that _D n (0) 2 + F n (0) 2 = D n (0) 2 + F n (0) 2 for all n < k 
and £> n (0) 2 + ^n(0) 2 = ^n- 2 (0) 2 + ^n- 2 (0) 2 for n > Jfe + 1. Hence £ P J0] = £p[0] + SE, 
where 5E is the energy of the pair of links forming the exceptional triple {4>k-i, 4>k, 0fc+i)- But 
5E = (4 + /i 2 )/4/i by Lemma 121 so 

E P [0]=Ep[0] + i^, (3.10) 
as was to be proved. □ 

If we choose = 0, the vacuum, in Proposition El then is a step-function kink located 
at k, of energy (4 + h 2 )/4h > 1. If is a kink, is a double kink whose energy exceeds 2, 
twice the energy of a single kink. If is an antikink, is a (step) kink-antikink pair. Clearly 
we can iterate the procedure to generate non-Bogomol'nyi solutions with arbitrary boundary 
values. All solutions so constructed have higher energy than the minimum within k+ (1 if 
k + — k_ is odd, otherwise), leading one to expect them to be unstable. Note also that the 
energy of any such solution diverges as h — > 0, as one would expect given that the continuum 
sine-Gordon model has no analogous static solutions. A static double kink, kink-antikink and 
triple kink are depicted in figure El 

In fact, one can prove that all non-Bogomol'nyi static solutions are unstable, independent 
of any particular construction. Crucial to the argument is the fact that every finite energy 
static solution has Bogomol'nyi "tails" asn^ ±oo. More precisely: 

Proposition 4 Let be a finite energy static solution. If has no exceptional triples, it is 
either a vacuum, a kink or an antikink. If has one or more exceptional triples, there exists 
N > such that D n = F n for all n > N or D n = —F n for all n > N, and D n = F n for 
all n < —N or D n = —F n for all n < —N. Hence coincides with a vacuum, a kink or an 
antikink for n > N, and likewise for n < —N. 

Proof: The proof rests on a pair of factorizations due to Ward jT^j. First, 
AdEp 1 



h dcbr, 2 



-COS-^n+l - <pn-l, 



1 

— sin -(0„+i - 20„ + n _i) - sin -(0 n +i + 20 n + n _i 

(311) 



S 




\ \ \ \ \ \ 



Figure 2: Some non-Bogomol'nyi static solutions of the spacing h* = 2/y/3 TDSG system 
(left), and their corresponding dipole multikinks (right): a double kink (a), a kink-antikink 
(b), and a triple kink (c). 



so every unexceptional triple in the static solution satisfies 



4 1 

^2 sin - (4>n+i - 20 n + 4>n- 



l) - Sin -(0 n +l + 2( Pn + <t>n-l) = 0. 



Second, 



(3.12) 



A( J D^ 1 -F n 2 _ i ; 



sm 



>n-\, 



>n+l 



20„ + 



>n-lj 



sm ■ 



>n+l 



+ 20 n + 



(3.13)" 

Now since Ep[(p] is finite, may have at most finitely many exceptional triples, by Lemma 
121 so there exists N > such that no triple {4> n -ii <f>n, <f>n+i) is exceptional for \n\ > N. But 
then for all n > N, ()3.12|) holds, so ()3.13|) implies that D\_ x — F^_ x — C, a constant. But 

and we conclude that = F% for all n > N. Similarly 
has no exceptional triples, the same argument shows that 



E„Pn + i ? n)<00, SO D n ,F n 



D 2 n = Fl for all n < -N. If 
D 2 n = Fl for all n. 

It remains to rule out the possibility (for n > N) that -D„_i 



while D r , 



or 



D r , 



-F n -x while D n = F n , for we can then conclude that <j) uniformly satisfies either ()3.7|) 



or (J3.8|) on its right tail. Assume the contrary. Then 



tan ■ 



tan ■ 



(3.14) 



so 4> n -i = n+ iinod27r. But this contradicts (|3.2jl . unless = 4> n+ i = or 7rmod27r. In 
this case, the right tail of is a vacuum since and tt are fixed points of both ()3.7j) and (|3.8)1 . 
Hence either satisfies ()3.7j) for all n > N, or it satisfies ()3.8j) for all n > N, or both (the 
vacua). The left tail, n < —N, is handled similarly. Again, if has no exceptional triples 
then either D„ = F„ for all n, or D„ = —F„ for all n, so is a kink, antikink or vacuum. □ 



9 



Theorem 5 Let be a finite energy static solution other than the kink, antikink or vacuum. 
Then is unstable, in the sense that there exists a6f 2 such that Hess^w, u) < 0. 

Proof: By Proposition HI has at least one and at most finitely many exceptional triples We 
may assume, without loss of generality, that the leftmost exceptional triple is (0_i,0o,0i), 
so that D 2 = F 2 for n < — 1. Let <p a,t be the 1-parameter family (parametrized by a) of 
variations (parametrized by t) of 0, defined by 



ia,t 
n 



J n n£{-l,0} 

l )- 1 + at n=-l (3.15) 
Jo + 1 n = 0, 



(note that a '° = for all a), and let w Q = d<p a,t / dt\ t= o £ We will show that there exists 
an a* such that 



Hess,;(//"\ //"•', 



E P [0 Q *'*] < 0. (3.16) 



t=o 



Note that w a = ae_i + eo, where {e n \ n G Z} is the standard basis for £ 2 , so 

Hess^w", m°) = a Q + a%a + a 2 a 2 =: p 2 (a), (3-17) 



E P [<p u > t \=0, (3.18) 



some quadratic polynomial in a. 

Now, 0°'* only varies by changing O , which affects only the energies of the links (0^, 0q *) 
and (0Q i ,0 1 ' i ). But these two links constitute an exceptional triple for all t, so the sum of 
their energies is (4 + h 2 )/4h, independent of t by Lemma 121 Hence Ep[<ft 0,t } = Ep[<j)} for all t, 
so 

d 2 

a = Hess^w ,?/ ) = — 

at t=0 

and p2{ot) = a\ct + 02a 2 . It thus suffices to show that p2(—l) 7^ £> 2 (1) ; since then a,\ 7^ so 
P2(oe) takes negative values close to 0. 

The variation <ft a,t changes only 0_i and 0o, so 

E P [(j) a,t \ = const - — {cos(0_i - 0_ 2 + at) + cos(0 o - 0_i + (1 - a)t) - cos(0 o - 0-i + £)} 

/i 

— - {cos(0_i + 0_ 2 + at) + cos(0 o + 0_l + (1 + at)t) - cos(0 o + 0-1 + t)} 
8 

1 h 
Pi 1 ) = 7TT {cos(0-i - 0_ 2 ) - cos(0 o - 0-l)} + - {cos(0_i + 0_ 2 ) + 3 cos(0 o + 0-l)} , 
2h 8 

1 ft, 

= — {COS(0_1 - 0_ 2 ) + 3 COS(0 O - 0-l)} + - {COS(0_1 + 0_ 2 ) ~ COs(0 O + 0-l)} 

2n 8 

^2a x = p(l)-p(-l) = -^cos(0o-0_i) + ^cos(0o + 0-i). (3.19) 

n 2 

But recall that D 2 _ 1 = F 2 1 , which implies that 



h 2 

cos(0 o - 0- X ) = 1 + — [cos(0 o + 0-i) - 1], (3.20) 
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and hence 

4 - h 2 

a x = —g- < 0. (3.21) 
Thus Hess^(w a , u a ) < for a sufficiently small and positive, so is unstable. □ 

We will now deduce from Theorem |3] that all non-Bogomol'nyi static solutions are linearly 
unstable. If the configuration space of our model were finite dimensional, this would follow 
immediately. However, in the infinite dimensional case it is necessary to proceed carefully. 
Let : £ 2 — ■» £ 2 be the Jacobi operator associated with the bilinear form Hess^, that is, the 
unique linear map such that 

(u, J+v)fl = Hess (w, v). (3.22) 

Then the linearization of ()2.6|) about the static solution <fi may be thought of as a second order 
linear ODE on £ 2 , namely, 

v + j* v = 0. (3.23) 

We think of v(t) G £ 2 as a perturbation of 0. If there exist perturbations which grow without 
bound (in £ 2 norm), we say that <fi is linearly unstable. Such a perturbation certainly exists 
if has a negative eigenvalue, since we can perturb in the direction of the corresponding 
eigenvector. It is not hard to show from Theorem El that spec contains at least one negative 
real number. However, we cannot conclude immediately that this number is an eigenvalue, 
since in infinite dimensions it may lie in the continuous or residual part of spec J^. A more 
detailed analysis of spec is required. We begin with two lemmas giving basic properties of 

Lemma 6 is bounded and self-adjoint. 

Proof: Clearly is symmetric by definition, so self-adjointness will follow from boundedness. 
As before, let {e n \ n 6 Z} be the usual basis for £ 2 and J% m := (e n , J^e m ). Then 

j, = 

" n.m. r\ i r\ i 



4 

— (cos(0 m - (j) m -i) + COs(0 m - 4> m+ l)) + COs(0 m + 4>rn-x) + COs(0 m + <j>m+l) 



4 

— cos(0 m - 4> n ) + cos(0 m + 



(3.24) 



Hence (J^v) n = c n v n + d n _iu n _i +d n+ iv n+ i, where c n , d n are sequences bounded independent 
of 0: \c n \, \d n \ < (4 + h 2 )/Ah. Thus 

ll^ll 2 = ^{CnVn + d n -lV n -l + d n+l V n+1 ) 2 

,2X2 'A + h 2 ^ 2 



^ 3 (^r) E^+^-i+^i)^^) in 



ii 



as was to be proved. □ 

It follows immediately from Lemma El that spec is real, bounded and consists of eigen- 
values and continuous spectrum only (there is no residual spectrum) [Tlj . 

Lemma 7 (Asymptotic positivity) Let <p be a finite energy static solution, and for each 
N G Z + ; let Vn = {v G I 2 \ v n = for all |n| > iV}. Then for all e > 0, there exists N such 
that 

(v,J*v)>(^-e)\\v\\ 2 

for all v G . 

Proof: By Proposition^ there exists Ni > such that coincides with a vacuum or (anti)kink 
located at b + for n > N\, and a vacuum or (anti)kink located at 6_ for n < —N\. Given the 
exponential decay properties of (anti)kinks ()2.7j) . (j3.9|) . it follows that 



COS(0 n+ i 



±<j> n ) -1| < Ce~ 2a|n| (3.25) 



for all \n\ > Ni + 1, where C > depends on b± and h, but not iVj.. (Recall that a = 
log((2 + h)/(2 — h)).) Let J° be the Jacobi operator of the vacuum, 

T n 4 + h 2 4 - h 2 , , , 

J% = -jj-v - -^(r+v + t-v), (3.26) 

where t± are the unit shift maps, (r±v) n = f n ±i, and L = — J°. 

Choose e > 0, and let iV be any integer exceeding both Ni and ^- log (^r). Then for all 
v _L Vn, v n — for |n| < iV, so we see from (|3.24j) and ()3.25|) that 

\(v,lv)\ < y: {^f ce ' 2aNv - + ^if\ v - v --A 

m<-N ' 

v—v f 4 + /i 2 __ 2aN 2 4 - /i 2 . .1 

m>iV v y 



4 + /i 2 ,. 2 4-/i 2 ^. . . 

Ml H 77— 2^ I'WWll f Ce 



2aV 



4/i " " 4/i 
< lce- 2aN \\v\\ 2 < e\ \v\ | 2 

At 

since 2|v m v m +i| < v 2 n + It follows, therefore, that 

(u,J*t;) = (u, J° v ) + (v,Lv) > (v, J°v) - \(v,Lv)\ 

/ t0\ 1 1 1 1 a ~ "' 1 1 1 1 2 / \ i i i 1 2 

> (u, J v) — e\\v\\ = — |M| — — (u,r+«) 

ft, 

> || 2 - e\\v || 2 

since r* = r + , (v,r + ) < \ \v\ \ \ \t+v\ \ by Cauchy-Schwarz, and r + is an isometry. □ 
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The key point is that J* is positive on a subspace of finite codimension (namely V^), from 
which we will see that it can have no negative continuous spectrum. Therefore, if J* has some 
negative spectrum, as it does when <fi is non-Bogomol'nyi, it must have at least one negative 
eigenvalue. We are now ready to state and prove the required spectral result. 

Theorem 8 Let <p be a finite energy static solution other than the kink, antikink or vacuum. 
Then has a negative eigenvalue of finite multiplicity. 

Proof: Since J (we will drop the superscript) is bounded and self adjoint, its spectrum is 
contained in the real interval [Qi,Q2\ where 

Qi := mf , Q 2 := sup , 

\\ v \\ v^O \\V\\ 

and, furthermore, contains Qi and Q2 jToj . By Theorem^ Q% < 0, so it remains to show that 
Qi is an eigenvalue of finite multiplicity. 

Choose e G (0, h/2). Then by Lemma[7[ there exists N > such that (v, Jv) > (| — e)||t> || 2 
for all v G V£. Let P : i 2 -»• i 2 be orthogonal projection onto V N , P' := 1 - P, T := P'JP' 
and A = P'JP + PJP' + PJP. Then J = T + A, both T and A are bounded and self-adjoint, 
and T is, by construction, positive. Given any bounded sequence u % in £ 2 , Au % is bounded, 
and takes values in a finite dimensional subspace of £ 2 (namely Vjv © P'J^Yn))-, so has a 
convergent subsequence by the Bolzano- Weierstrass theorem. Thus A is compact, and hence 
trivially T-relatively compact [TBJ pl94]. Thus spec e J, the essential spectrum of J = T + A 
coincides with the essential spectrum of T [TBI p244], which is bounded below by | — e [T5] . 
Hence Qi is not in spec e J. But J is self-adjoint, so spec J\spec e J contains only eigenvalues 
of finite multiplicity fUJ p518]. Hence Q\ is a negative eigenvalue of finite multiplicity, as was 
to be proved. □ 

We note in passing that, since e could be chosen freely in the above argument, we get the 
extra information that the continuous spectrum of J is bounded below by \. 

Corollary 9 (Linear instability) Let <ft be a finite energy static solution other than the kink, 
antikink or vacuum. Then the linearized TDSG flow about (f> supports solutions which grow 
unbounded exponentially fast. 

Proof: Let u G i 2 be the eigenvector associated with the negative eigenvalue Q\ =: — v 2 . Then 
(I3.23J1 supports the solution 

v(t) = uexp(J^ut), (3.27) 
which grows unbounded exponentially fast. □ 

We close this section by noting that there is a two parameter family of static solutions for 
which every triple is exceptional: 

= mod 4 

= l m ° d A (3-28) 
= 2mod4 v ; 

= 3 mod 4. 
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These stand at the opposite extreme from the Bogomol'nyi solutions. Clearly they have infinite 
energy. 



4 Exact propagating kinks 

An important question is whether there exist moving kink solutions that propagate with 
constant velocity. In other words, in this section we are going to look for solutions that satisfy 
the following condition: 

4>n(t) = <j>{n — st) = cj)(z) , z = n — st . (4.1) 

Substituting this ansatz into (|2.fij) . we see that the profile <fi satisfies a nonlinear advance-delay 
ODE, namely 

4 - h 2 

s 2 (p"(z) = ^ 2 cos 4>{z) [sin cf){z + 1) + sin cf){z — 1)] 
4 + h 2 

sin (j)(z) [cos 4>(z + 1) + cos <j>(z - 1)] . (4.2) 



Ah 2 

It is very hard to prove rigorous results on the solutions of such ODEs, and comparatively 
little is known in general. Friesecke and Wattis jlZj have proved the existence of propagating 
pulses in FPU chains with superquadratic intersite potential. Somewhat closer to the current 
situation, Iooss and Kirchgassner [TH] have proved the existence, in linearly coupled oscilla- 
tor chains, of small amplitude travelling pulses with small oscillatory tails. Neither method 
developed in these papers applies directly to the travelling kink problem we seek to address. 

We will approach the problem numerically by seeking kink solutions of ()4.2j) using the 
pseudo-spectral method, a highly effective tool for finding travelling-wave solutions in lattices 
which has been developed in a number of papers fHJ 1201 121] The idea is to write as a basic 
kink (for example 2tan _1 e az ) plus a small correction 54>(z), write down a truncated Fourier 
series for S(p(z) and hence reduce ()4.2j) to a large algebraic system for the unknown Fourier 
coefficients, which can be solved using the Newton method. For technical details see |20[ I21j. 
We stress that this method is based not on lattice simulations but on iterative techniques and 
it enables us to find travelling-wave solutions with any given precision. 

Straightforward application of the pseudo-spectral method to the TDSG system shows 
that, as for conventional discretizations, moving kink solutions with oscillating tails, or so- 
called nanopterons, appear. Two such nanopterons are depicted in figure 0J The amplitude 
A of the oscillating tail depends on speed s and lattice spacing h. If A(s, h) = 0, one has 
a genuine travelling kink, since <ft satisfies kink boundary conditions. Such solutions usually 
have codimension 1 in the s, h parameter space, so they occur at isolated values of s if h is 
fixed, and vice versa. A heuristic explanation of this was provided by Aigner et al [221 • They 
think of the travelling kink profile as a heteroclinic orbit from to 7r and informally compute 
the dimensions of the stable and unstable manifolds of these fixed points by considering the 
linearization of (|4.2|) about them. Every independent oscillatory solution of the linearization 
cuts the dimension of each of these manifolds by 1, and increases the codimension of the space 
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of travelling kinks by 1. Applying their argument to the TDSG system leads one to expect 
that the codimension equals the number of non-negative solutions k of the equation 

s 2 fc 2 = l + i^sin 2 |, (4.3) 

obtained by demanding that <p = cos kz be a solution of the linearization of ()4.2|) . For fixed 
h, this number is 1 for s large, but grows without bound as s — > 0, suggesting that travelling 
kinks, if they exist at all, should disappear in the low speed limit. On the other hand, the s — > 
limit of the TDSG system is clearly very wn-generic, since the system supports a continuous 
orbit of static kinks, so it remains possible that the arguments of are misleading in this 
case. 




Figure 3: Profiles of the nanopteron solution u(z) of ()4.2|) for s = 0.2 and s = 0.5 (with larger 
oscillating tail). The inset shows details of the tail of the solution for s = 0.2. In both cases 
h = 1.3. 

To resolve this issue, we have systematically computed the nanopteron tail amplitude A as 
a function of s for various values of h, see figure EJ In no case do we see zeros of A, isolated or 
otherwise, at large s, so we actually find no travelling kinks in the portion of parameter space 
where they are expected to arise with codimension 1. The small s regime is more interesting. 
As s — > 0, the amplitude drops monotonically to a value so small as to be to within numerical 
tolerance. Certainly A(0) = since there is no PN barrier, and <fi(z) = 2 tan -1 e az solves ()4.2j) 
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at s = 0. This small s behaviour is very different from that of the conventional discrete 
sine-Gordon system (DSG), 



Wl-20n + ^_ s . n ^ (44) 



h 2 

included in figure El for comparison: A(s) clearly remains bounded away from for this system. 
For the TDSG system, it is possible that A(s) really does attain the value on some narrow 
interval [0, s*(/i)], so that exact travelling kinks exist for all speeds not exceeding s*(/i). More 
likely A(s) rises immediately from 0. This would be consistent with some formal asymptotics 
of Oxtoby et al computed in the context of exceptional discretizations of the 4 system |23j . 
However, the question is numerically inaccessible. All we can say is that the tail amplitude 
goes rapidly to zero as s — > 0, and that this is reflected in a high degree of practical kink 
mobility, as demonstrated in previous studies 012111], and by our own simulations, described 
below. 




s 



Figure 4: Amplitude of the oscillating tail as a function of the kink velocity s for h = 1 (o), 
h = 1.3 (V), h = 1.5 (©) and DSG (x) [given by Eq. (|4.4|) with h = 1]. Connecting lines are 
used as guides to the eye. 

In order to check the dynamics of the TDSG and DSG systems, we performed numerical 
simulations of (|2.6|) and (J4.4)) . The initial conditions were chosen by assuming translational 
invariance of the solution. In the case of the TDSG, we substitute b = st in (|2.7|) and compute 
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n (O),0 n (O): 

n (O) = 2tan- 1 e a(n - no) , (4.5) 

0n(O) = -r;(f>n(t)\t=o = -as/ cosh[a(n - n )] . (4.6) 

In the case of the conventional DSG, due to the absence of an explicit stationary solution, we 
take its continuum approximation n (t) = 4 tan -1 e h ( n - n o- st ) anc i compute n (O),0 n (O). Both 
equations have been simulated using the 4th order Runge-Kutta method on a finite lattice of 
size N. The position of the kink's centre of mass, defined as 

XM = , (4.7) 

2(4> N - (pi) 

has been plotted function of time in figure E3 




Figure 5: Temporal evolution of the kink's centre of mass for h = 1.3 for two discretizations of 
the sine-Gordon equation. Lines 1, 2, 3 and 4 correspond to the TDSG with initial velocities 
s = 0.001, s = 0.01, s = 0.1 and s = 0.5, respectively. Line 5 corresponds to the kink of the 
conventional DSG kicked with initial velocity s — 0.1. In the simulations, dissipation has been 
introduced at both ends of the chain. The left inset shows details of lines 3 and 4, and the 
right inset shows details of line 5. 

One can clearly see from this figure the difference between the two discretizations. In the 
TDSG case the kink moves freely along the lattice even when kicked with very small initial 
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velocities (lines 1, 2, and 3). By inspecting this figure it is easy to see that, for small initial 
velocities, the actual velocity of propagation almost coincides with the initial velocity, so that 
the kink suffers practically no radiative deceleration. The figure also shows the existence of 
a velocity threshold above which radiative deceleration becomes important: the kink given 
initial velocity s = 0.5 quickly slows down to approximately s — 0.1 (line 4). By contrast, line 
5 corresponds to the kink of the conventional DSG, kicked with an initial velocity s — 0.1. The 
kink propagates several sites and then finds itself trapped by the lattice. The same happens 
for smaller values of the initial velocity s. Thus, in the dynamical picture, the absence of the 
PN barrier for the TDSG manifests itself in the possibility, for all practical purposes, of free 
propagation of kinks with an arbitrarily small velocity. 

5 Conclusion 

We have shown that there is a concrete physical system modelled by an exceptional discretiza- 
tion of the sine-Gordon model, namely a uniform chain of electric dipoles, modelled by the 
TDSG system. Previous work [2] shows that such a chain admits a continuous translation 
orbit of static kinks, with no PN barrier, and that these kinks are highly mobile and exhibit 
strongly continuum-like behaviour (no trapping, low radiative deceleration etc.), despite being 
inherently very highly discrete. We have constructed static multikink solutions of the system 
and proved that all finite energy non-Bogomol'nyi static solutions (i.e. all except the kink, 
antikink and vacua) are unstable. 

A key step in the analysis of static solutions was the observation that, whilever the solution 
has no exceptional triple, it "conserves" the two-point function — F% as n "evolves" along 
the lattice. This is strongly reminiscent of Kevrekidis's method so one might ask whether 
the TDSG system conserves both energy and momentum (it falls outside the class of models 
considered in jTj). The answer appears, unfortunately, to be no. From (|2.6|) . (j3.11|) and ()3.13|) 
it follows that 



for any finite energy solution, but there is no obvious way to write this as a total time derivative. 

We have investigated the possibility that the TDSG system admits exact constant-speed 
propagating kinks, finding only nanopterons, that is, kinks with spatially oscillatory tails. 
This is in contrast to the results of Oxtoby et al on exceptional discrete 4 systems, who 
find isolated velocities for which the tail amplitude vanishes exactly so that exact propagating 
kinks exist |2H|. In the TDSG case, the tail amplitude vanishes rapidly in the zero velocity 
limit, and may, in fact, remain exactly zero on a small velocity interval, but the resolution 
of this issue is beyond the capability of our numerics. If so, the travelling kinks would form 
a codimension family in a regime where, if they exist at all, they should be isolated, in 
contradiction to the heuristic barrier of Aigner et al [22j. This is perhaps unlikely [23] , but 
should not be discounted altogether given the other highly ungeneric features of the system. 
Whether or not exact propagating solutions exist, numerical simulations confirm that the kink 
can move freely with arbitrarily small velocity, so that kinks enjoy a high degree of practical 
mobility, in marked contrast with the conventional DSG. It is possible that exact travelling 




(5.1) 



n 
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kinks with a structure more complicated than that accommodated by the ansatz (|4.1|) exist. 
This problem requires a separate investigation. 

To our knowledge, this is the first known example of an exceptional discrete system which 
models a genuine physical system. One could object that the model takes into account only 
the interactions between nearest neighbour dipoles. This is true, although this is a standard 
approximation made when modelling lattice systems. Including longer range inter-dipole forces 
will presumably destroy the exact continuous translation orbit of static kinks and introduce 
a small PN barrier. Since longer range forces are much weaker (note the |r| -3 in (j2.1|0 one 
would hope that the PN barrier is so small as to be negligible in practice. To estimate how 
small, we have added to Ep the energy due to next-to-nearest neighbour dipole pairs, 

E 'p = I 4 J2 { jj? sin2 - (pn-i) + sin 2 ^(0 n +i + <f>n-i) \, h := Vl2, (5.2) 

in natural units, and numerically minimized Ep + E' P for site centred and link centred alter- 
nating kinks using a gradient flow method. We find that site centred kinks now have very 
slightly lower energy. The PN barrier 

E(b=\)-E(b = 0) 

E{b^0) < °-° 09% - (5 ' 3) 

This is absolutely tiny in comparison with conventional discrete systems, and strongly suggests 
that the inclusion of only nearest neighbour interactions in the model is a sensible approxi- 
mation. 
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